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Foliage profiles retrieved from a scanning, terrestrial, near-infrared (1064 nm), full-waveform lidar, the Echidna 
Validation Instrument (EVI), agree well with those obtained from an airborne, near-infrared, full-waveform, 
large footprint lidar, the Lidar Vegetation Imaging Sensor (LVIS). We conducted trials at 5 plots within a conifer 
stand at Sierra National Forest in August, 2008. Foliage profiles retrieved from these two lidar systems are closely 
correlated (e.g., r = 0.987 at 100 m horizontal distances) at large spatial coverage while they differ significantly 
at small spatial coverage, indicating the apparent scanning perspective effect on foliage profile retrievals. Also we 
noted the obvious effects of local topography on foliage profile retrievals, particularly on the topmost height re- 
trievals. With a fine spatial resolution and a small beam size, terrestrial lidar systems complement the strengths 
of the airborne lidars by making a detailed characterization of the crowns from a small field site, and thereby 
serving as a validation tool and providing localized tuning information for future airborne and spaceborne 
lidar missions. 

Published by Elsevier Inc. 


1. Introduction 

The vertical distribution of foliage area per unit volume above the 
ground, referred to as the foliage profile, shows great promise and ca- 
pacity for the estimation of aboveground biomass (Drake et al., 2002; 
Lefsky et al., 1999; Stark et al., 2012), for the distribution of transmitted 
directed light inside the canopy (Parker et al., 2001 ), and for the estima- 
tion of gross primary productivity (GPP) and net primary productivity 
(NPP) (Kotchenova et al., 2004; Toda et al., 2009). However, direct 
ground measurements of the vertical foliage distribution, such as with 
stratified clipping of biomass samples, are generally destructive, time 
consuming and labor intensive. Similarly, indirect techniques such as 
the point-quadrat method (Wilson, 1960, 1963, 1965) are also labori- 
ous and time consuming, especially at the forest stand scale. As a result, 
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considerable research has been devoted to measuring foliage profiles 
with remote sensing approaches that are particularly suited for 
large-area sampling and mapping of spatially explicit estimates of fo- 
liage profiles (Harding et al., 2001; Jupp et al., 2009; Lefsky et al., 
1999; Strahler et al., 2008; Tang et al., 2012; Zhao et al., 2011). 

Most airborne and spaceborne optical remote sensing systems pro- 
vide estimates of the horizontal distribution of canopies (e.g., Landsat, 
MODIS), but have difficulty in characterizing the vertical distribution 
of canopy elements (Harding et al., 2001; Sun et al., 2008). Only light 
detection and ranging (lidar) systems record the distance between 
the sensor and a target, and the resulting distribution of backscattered 
energy over height directly reflects the three dimensional vertical struc- 
ture of vegetation elements. Therefore, lidar scanning systems are being 
widely used for foliage profile retrievals in recent studies (Harding et al., 
2001 ; Hilker et al., 2010; Jupp et al., 2009; Lefsky et al., 1999; Strahler et 
al., 2008; Zhao et al., 2011). Foliage profile mapping over a large area 
generally requires the deployment of a discrete-return lidar scanning 
system or a full-waveform lidar scanning system on an airborne or 




F. Zhao et al / Remote Sensing of Environment 136 (2013) 330-341 


331 


spaceborne platform (e.g., the airborne full-waveform lidar, Laser 
Vegetation Imaging Sensor (LVIS), airborne discrete-return lidar, 
TopEye and the spaceborne full-waveform GLAS system aboard on 
the Ice, Cloud, and Elevation Satellite (ICESat)). Compared with a 
discrete-return lidar system, a full-waveform lidar scanning system 
shows much greater potential for deriving information on the details 
of the foliage elements over height. First of all, full-waveform lidar scan- 
ning systems typically have much larger footprints (e.g. 10-60 m) that 
are commensurate with the scale of structure variation in the canopy 
(Harding et al., 2001). Second, the footprints of discrete-return lidar 
scanning systems are generally kept small (e.g., 10-30 cm) to focus 
lidar energy for higher signal-to-noise ratios in backscattered energy. 
As a result, depending on the sampling density of the lidar pulses, 
discrete-return systems may fail to characterize the lower canopy in a 
relative dense forest stand, or miss some tree tops. 

It must be noted, however, that the ability of large-footprint, 
full-waveform lidar to retrieve foliage profile over sloped regions is 
limited, due to the complex interactions of lidar signals from vegeta- 
tion and ground surface (Yang et al., 201 1 ). Local topography affects 
the foliage profile retrieval directly by stretching waveforms, and the 
stretching effect magnifies with footprint size, slope, and off-nadir 
pointing angle (Lee et al., 2011). For example, within 20° of slope, 
given a height of roughly 25 m and a footprint size of 25 m, the er- 
rors for vegetation height range from — 2 m to greater than 20 m 
using an extended geometric optical and radiative transfer model 
(Lee et al., 201 1 ), implying that the effect of slope on foliage profiles 
needs to be carefully evaluated when comparing different lidar sys- 
tems. Meanwhile, Ni-Meister et al. (2001) and Kotchenova et al. 
(2004) reported that the retrieved foliage profiles from airborne 
lidar scanning systems differ from the actual foliage profiles by the 
mean projection of normals on the direction of laser beam because 
the airborne system can only provide near nadir-view sampling. 
Therefore, foliage profiles retrieved from airborne or spaceborne 
platforms are generally termed as “apparent foliage profile” or “can- 
opy height profile” while “actual foliage profile” or “foliage profile” is 
used for foliage profiles derived from lidar system sampling over the 
entire zenith region ranging from 0 to 90°. 

Ground-based, full-waveform lidar scanning systems (e.g., the 
full-waveform lidar scanning instrument Echidna Validation Instru- 
ment (EVI)) capture both the upper hemispherical view of the canopy 
and the detailed description of local topography, and these systems 
have been used successfully to retrieve actual foliage profile retrievals 
(Hosoi and Omasa, 2007; Jupp et al., 2009; Strahler et al., 2008; Yao et 
al., 2011; Zhao et al., 2011). Such detailed characterization of forest 
structure from terrestrial lidar systems should provide more accurate 
ground-truth foliage vertical structure at the plot scale, and compar- 
ing foliage profiles retrieved from a terrestrial full-waveform lidar sys- 
tem with those of an airborne full-waveform lidar system would help 
understand the mechanism leading to the difference between apparent 
foliage profile and actual foliage profile, and also the effect of local to- 
pography on foliage profile retrievals. In addition, in light of the wide 
use of downward-scanning, full-waveform lidar systems for foliage pro- 
file retrievals, for example the airborne full-waveform Laser Vegetation 
Imaging Sensor (LVIS) and the spaceborne full-waveform GLAS system 
aboard the Ice, Cloud, and Elevation Satellite (ICESat), such a compari- 
son of foliage profile retrievals is of particular interest. 

The scanning perspective is a crucial issue in comparison between 
terrestrial lidar foliage profiles and airborne lidar foliage profiles pri- 
marily because terrestrial lidar systems feature a radial perspective of 
observations within a limited range, and thus their capability to charac- 
terize the forest is heavily influenced by the forest density as well as the 
sampling design used. Fig. la, b shows a top-down view of a point cloud 
for the center scan and for five merged scans for a red fir stand sampled 
in the Sierra National Forest. These two figures demonstrate a pattern 
where points are densely distributed around the center, and this density 
decreases gradually away from the center. Fig. lc and d shows the point 


density distribution histograms corresponding to the center scan 
and to the five merged scans, corroborating the pattern shown in 
Fig. la, b. In contrast, the airborne lidar scanning systems are generally 
nadir-looking, and sample the forests in an equally spaced periodic 
manner. Therefore, the impact of the varying scanning perspectives on 
the retrieved foliage profiles must be considered when comparing re- 
sults from upward- and downward-looking lidar systems. 

In this study, we compare foliage profiles from a downward-looking, 
airborne, full- wave form lidar scanning system (i.e., LVIS) with those 
from an upward-looking, terrestrial, full-waveform lidar scanning sys- 
tem (EVI) for the foliage profiles. Specifically, we focus on 1) assessing 
how the foliage profile retrievals are impacted by the varying scanning 
perspectives of EVI and LVIS systems; 2) evaluating the effect of surface 
topography on the foliage profile retrievals from EVI lidar system; and 
3) finally comparing the foliage profiles retrieved from the EVI and 
LVIS systems. 

2. Methodology 

2.1. Study area 

Data from the EVI and from ground measurements were acquired 
from a high-elevation red fir stand in the Sierra National Forest, CA, in 
July, 2008. The sample area covers 1 ha in a square of 100 by 100 m 
and is subdivided into nine plots of 33 x 33 m 2 . The sample layout 
was oriented in azimuth so that plot sides followed the primary direc- 
tions downslope and across slope. We measured tree characteristics 
and leaf area index (LAI) with both the LAI-2000 and with hemispheri- 
cal photos (13 observations per plot). Five EVI scans were acquired- 
one in the center of the center plot and four at the center of the sur- 
rounding corners of that plot (Fig. 2a). The black cross signs represent 
the scan positions of EVI instrument. 

LVIS data were collected during September, 2008, within a few 
weeks of the EVI scans. LVIS Ground Elevation (LGE), LVIS Canopy Top 
Elevation (LCE), and LVIS Geolocated Waveform (LGW) data were 
downloaded from the Goddard Space Flight Center's LVIS website 
(lvis.gsfc.nasa.gov) (Blair et al., 2006). The “zg” (i.e., the mean elevation 
of the lowest mode within the waveform) field in the LGE file was used 
as the ground elevation for each waveform. Note that both horizontal 
accuracy and vertical accuracy for the sub-canopy topography were 
found to be approximately 2 m (Ict) for a dense, 98 to 99% closed trop- 
ical forest in Costa Rica (Blair et al., 1999). 

To match the center of EVI footprints with those of LVIS over-flights 
in August 2008, the geographic coordinates of each EVI scanning point 
were recorded using a hand-hold Garmin GPSMap60 CSx with accuracy 
of ± 10 m and then used to identify LVIS waveforms falling into the EVI 
ground footprints. In each ground plot, all live trees and snags 10 cm or 
greater in diameter at breast height (DBH) were targeted for structural 
measurements, including species, DBH, crown diameter, tree height, 
and height to the bottom of the partial crown. Fig. 2b shows a sample lay- 
out of the LVIS shots with the EVI ground point cloud as the background. 
The color gradient reflects the vertical height gradient of these forest 
plots. Note in Fig. 2b that the EVI scan center may not precisely match 
a single waveform, due to the inherent systematic errors of both the 
GPS device mounted on the LVIS system and the GPS device used to col- 
lect geographic coordinates of each EVI scanning point. (Fig. 3) provides 
a side view of a point cloud before slope correction and its corresponding 
view after slope correction. 

2.2. LVIS data processing 

LVIS is a full-waveform-digitizing system and records the vertical 
distribution of nadir-intercepted surfaces at a 30 cm vertical resolu- 
tion. The LVIS instrument flew aboard the aircraft at 7 km above 
ground level, and recorded footprints illuminated within a =t 7° field 
of view. For the 2008 flights, these footprints had a nominal diameter 
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Fig. 1. Scanning perspectives of EVI data, a) A top-down view of a red fir site in the Sierra National Forest center scan; green tones from dark to light indicate heights from low to 
high, b) A top-down view of five merged scans, c) A histogram of point density over range for the site 305 center scan, d) A histogram of point density for the five merged scans of 
site 305. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


of 20 m, and therefore the canopy vertical structures should be well 
characterized with sufficient LVIS shots (Sun et al., 2008). The actual 
footprint size varies with the altitude of the aircraft, the incidence 
angle of the pulse, and the elevation of the target. 

The methodologies for obtaining canopy height profile using air- 
borne waveform-recording laser altimeters require a series of as- 
sumptions (Harding et al., 2001; Lefsky, 1997). First, the horizontal 
distribution of the canopy elements is assumed to be random with a 
constant leaf inclination distribution. Second, the ratio of woody to 
leaf surface area and the woody and leaf reflectance is constant as a 
function of height. Finally, multiple scattering contributes only slightly 


to the backscattered waveform signals. The basic steps for estimating 
canopy height foliage profile originated with MacArthur and Horn's 
(1969) method for a reflex camera with the following expression: 


LAI, 


In (\/\ 

h 2 -hi 


( 1 ) 


where N hi and N hi are the number of points on which the lowest 
leaves are higher than hi and h 2 , respectively, and LA/ h2 _ hi is the leaf 
area density between and h 2 . Given the number of points from 



Fig. 2. a) Layout for EVI scans at Sierra National Forest; b) layout of LVIS shot center (black dots) laid upon virtual EVI center scan region of 20 m radius (red circle) with the down- 
ward view of EVI ground point cloud (dark green) at site 305 in the background. (For interpretation of the references to color in this figure legend, the reader is referred to the web 
version of this article.) 
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Fig. 3. EVI lidar point cloud for site 305. a) Before topographic correction; b) after to- 
pographic correction. 

the ground to the top of the canopy, the foliage area volume density 
can be obtained by dividing the total LAI by proportion of each layer 
to the total. The waveform-recording laser instrument, however, mea- 
sures the backscattered energy distributed over height, rather than 
providing the number of intercepted points, and therefore the lidar 
method needs adjustments for foliage area retrievals. The adjustments 
include smoothing the signal, identifying the background noise level, 
differentiating the ground and canopy returns, adjusting the return 
amplitude to account for difference in reflectance, computing a height 
distribution of canopy closure and transforming the measure to a cu- 
mulative distribution of plant material (Harding et al., 2001 ). 

However, as mentioned, the resulting foliage profile is not exactly the 
actual foliage profile because it represents solely the nadir-intercepted 
surfaces. Ni-Meister et al. (2001) reported that the retrieved canopy 
height profile deviates from the actual foliage profile by the mean projec- 
tion factor expressed in the following expression. 


FHP retr (z) = FHP act (z)G(z,0 L ) (2) 

where G(z, 0 L ), hereafter referred to as the G function, is the mean 
projection of leaf normals on the direction of laser beam 0 L at height 
z, and FHP retr (z ) and FHP act ( z) are the retrieved and actual foliage 
area volume densities respectively. The G function varies with zenith 
angle, and its variation depends on the angle distribution of the facets 
among the scattering elements (Jupp et al., 2009). The G function is 
only needed for leaf area index and foliage profile retrievals when 
limited zeniths are sampled. When all zeniths are sampled, LAI and 
foliage profiles can be retrieved without specific knowledge of the G 


function (Jupp et al., 2009). Kotchenova et al. (2004) chose to use a 
field-measured canopy area index (CAI) to normalize the retrieved 
foliage area profile. 


FHP(z) = ( CHP retr (z)/CAI retr )LAI (3) 

where CHP retr is the retrieved canopy height profile, CAI retr is the re- 
trieved canopy height index and LAI is the field measured LAI. 

In this study, raw LVIS waveforms were processed following the pro- 
cedure described in Harding et al. (2001 ), and then the retrieved CHPs 
were averaged at varying spatial coverages, followed by normalization 
using Kotchenova's approach. All LVIS waveforms within a specified 
range from the EVI scanning center were extracted and added together 
by aligning them with waveform ground peaks, resulting in mean fo- 
liage profiles. 

For the LVIS system, the height of the canopy height profile (CHP) 
derived from a single waveform is referenced to the start of the ground 
peak, and as a result, the between-footprint topography effect is mini- 
mized when several foliage profiles are averaged over an area for a 
mean foliage profile. However, the within-footprint topography effect 
remains inherent in each waveform, and cannot be removed directly. 
As a result, in this study, the foliage profiles retrieved from LVIS wave- 
forms are corrected for the between-footprint effect, but still retain 
the within-footprint topography effect. Fig. 4 shows a DEM for site 
305 interpolated from the mean elevation of the lowest mode within 
the waveforms recorded in the LGE file. 

Locational errors associated with both the recorded geographic co- 
ordinates of the EVI scanning centers and those of LVIS waveforms 
using GPS devices contribute to a slight mismatch between single LVIS 
waveforms and EVI scanning centers. Therefore, the production of a 
mean foliage profile over a larger region reduces the impact of the 
mismatch of foliage profiles derived from upward-looking EVI and 
downward-looking LVIS. 

To match the foliage profile from the upward-looking EVI with 
LVIS, we produced mean foliage profiles within circular regions 
with varying horizontal ranges (radii) from the center point, ranging 
from 20 to 100 m. Hilker et al. (2010) showed that in a restricted 
geographic coverage (i.e., 30 x 30 m), foliage profiles from airborne 
discrete-return lidar match well with those from EVI. Here, we chose 
a circular-shaped pattern of areas of increasing horizontal range that 
better matches the spherical scanning of the EVI and the circular 
footprint of the airborne LVIS scanner. 

2.3. Foliage profiles retrieved from EVI full-waveform data 

For the EVI system, if the probability of canopy gap (hereafter re- 
ferred to as Pgap) at the top of the canopy for a range of zenith angle 
rings is known, the LAI and vertical foliage profile can be readily 
obtained (Zhao et al., 2011). The mean profiles, excluding the zero ze- 
nith (or first ring) since it is the most variable, are averaged using the 
solid angles subtended by the rings as weighting. 

Previous studies of the topographic effect on LAI estimated from ver- 
tically upward pointing hemispherical photographs have shown that 
slope increases the gap fraction downslope while it decreases gap frac- 
tion upslope. Walter and Torquebiau (2000) proposed replacing zenith 
angle 6 with incidence angle y (relative to the normal to the slope), 
which is defined as: 

cosy = cos0 cos/3 + sin0 sin/3 cos (<p— <p) (4) 

where 6 is the zenith angle (relative to the flat ground), /3 is the ground 
slope, <p is the azimuth of incidence angle and <p' is the azimuth of 
ground slope. In this condition, the path length is no longer held con- 
stant for all azimuth directions. 
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Fig. 4. DEM derived from LVIS waveforms. Green color from dark to light indicates height from low to high. (For interpretation of the references to color in this figure legend, the 
reader is referred to the web version of this article.) 



To remove the local topographic effect on foliage profiles retrieved 
from the EVI full-waveform data, it is essential to identify how the 
ground returns fit along a regular slope. As these ground shots have 
somewhat unique characteristics (e.g. single returns, strong peak signa- 
tures, and primarily distributed in the lower hemisphere), they can be 
easily identified for further slope fitting. Once the single ground shots 
are identified with a form of point cloud, the slope can be fitted mathe- 
matically in 3-D space. 

We fitted a least-squares plane to the elevation of the ground shots 
using an orthogonal distance regression that minimizes the sum of the 
distances of all points to the fitted plane. The slope and azimuth can 
then be derived based on the fitted equation. 


Ax + By + Cz + D = 0 (5) 

where A, B and C are the components of the plane. Once the slope 
angle and azimuth angle were determined, we used the methodology 
proposed by Walter and Torquebiau (2000) to see how the foliage 
profile estimates were affected by local topography. 

Compared with the MacArthur and Horn's (1969) method used to 
retrieve the foliage profiles, the gap fraction-based foliage retrieval pro- 
files can only be processed at the scanning point level. To match the spa- 
tial scale, we can compare foliage profiles of a single center scan with 
those profiles obtained within a small spatial coverage area (e.g., 30 m 
in radius), while for larger spatial coverages, we need to use the mean 
foliage profile retrieved from five EVI scans to match the mean foliage 
profiles from LVIS within the same coverage. 

2.4. Foliage profiles retrieved from EVI point cloud data 

As opposed to airborne lidar scanning systems, with a fine spatial 
resolution and small beam size, terrestrial lidar systems can comple- 
ment the weaknesses of airborne lidars by making a detailed character- 
ization of the crowns in a small field site, making more accurate foliage 
profiles possible. For example, approximately 400,000 lidar pulses will 
exit the canopy over an area of 0.85 ha (assuming a 30-m canopy 
height) at varied view angles for EVI lidar scanning system. With 5 
scans and about 2 million shots, nearly all the foliage region in 1-ha 


site will be well sampled. Also, more than 1 00,000 shots are directed 
to the ground from each scan, and as a result returned pulses from 
five merged scan can provide a quite detailed characterization of local 
topography (Fig. 5). 

Although gap fractions can be obtained using a gap probability-based 
method at varying zenith regions or a gap fraction method over height 
(Hilker et al., 2010; Hosoi & Omasa, 2006), the foliage profile can also sim- 
ply be retrieved from a point cloud using the following two expressions 
(Jupp et al., 2009). 


a(- log (p(z))) 

m = A — 11 ( 6 ) 

where P gap {z) is the gap fraction measurement at height z; and /(z) 
represents the foliage area volume density at height z. In this study, 
the gap fraction can be estimated as the total number of lidar hits 
up to a height z relative to the total number of lidar shots emitted 
out of the EVI instrument (Coops et al., 2007). 



where #Z is the number of hits from the default EVI height 1.7-m up 
to a height z above the ground; N is the total number of shots emitted 
up to the sky. To avoid any hits that are intercepted by upslope local 
topography, only shots with zenith angles less than 85° were chosen. 

Although the EVI can record hits beyond a 1 00-m range, most of the 
lidar shots are intercepted in a closer limited range, as shown in Fig. lc, 
d. This is particularly true of the shots reflected back from the ground 
because the EVI instrument is mounted on a tripod at approximately 
1.7 to 2.0 m. As a result, the best characterization of local topography 
is achieved by merging several EVI scans together. 

To correct for the topographic effect on foliage profile retrievals from 
the EVI instrument, the ground points were first identified, and then 
fitted to a three-dimensional surface. The heights of other points were 
then adjusted accordingly using this 3-D surface. 
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The EVI point cloud used in this study is transformed directly from 
the EVI full-waveform data by identifying the waveform peaks within 
each single waveform as hits; each of the hits is then recorded with its 
geographic coordinates, which provides much more flexibility in re- 
trieving foliage profiles with varying spatial coverages. To evaluate 
how the difference in scanning perspective impacts the foliage profile 
comparison, this study chose to provide the foliage profiles from vary- 
ing circular regions, in the same way we retrieve mean foliage profile 
from LVIS. By investigating the foliage profile comparison at increasing 
circular ranges from a center point, a better understanding of the effect 
of the scanning perspective on the foliage profiles from these two 
sources is achieved. 

3. Results 

3.1. Foliage profiles retrieved from LVIS 

Fig. 6 shows the LVIS-derived apparent foliage profiles as observed 
with circular ranges varying from 20 to 100 m for a single red fir site 
305. Generally, the foliage profiles retrieved from closer horizontal 
ranges have relatively smaller FAVD values for the lower canopy (i.e. 
12-18 m) and upper canopy (40-60 m), compared with those re- 
trieved at longer ranges, while for the middle canopy (i.e. 30-40 m), 
they display a relatively higher FAVD. This variation demonstrates 
both the horizontal and vertical heterogeneity of the foliage distribution 
in this stand. Horizontally, as more single foliage profiles were aver- 
aged, the shape of mean foliage profiles became smoother, and this re- 
flects a comprehensive view of this forest stand. Vertically, the canopy 
top heights of the foliage profiles retrieved at close ranges are higher 
than those of the foliage profiles at farther ranges, reflecting the 
mixed effect of slope and vertical heterogeneity over larger spatial re- 
gions. The difference in top height between the 20 m and 100 m ranges 
approaches 10 m, indicating the necessity of addressing the local topo- 
graphic effects on the foliage profile retrievals, and also the need to 
match the foliage profiles spatially in a rigorous manner. Yang et al. 
(2011 ) assessed the generalized impacts of surface topography on veg- 
etation lidar waveforms using an extended geometric optical and 

a 3 (k 
20 ' 



-100 -100 

Y (m) x < ra > 


Fig. 5. Digital elevation model (DEM) derived from EVI point cloud, a) DEM derived 
from the EVI ground returns; b) fitted to a plane slope surface. 



Fig. 6. LVIS-derived foliage profiles with varying spatial coverages, shown by horizon- 
tal radii (20-100 m) of coverage areas. 

radiative transfer model. However, to remove the slope effect inherent 
in a single waveform obviously requires a more detailed description of 
forest structure and local topography. 

3.2. Foliage profiles retrieved from full-waveform EVI observations 

Fig. 7a, b illustrates the specific local topographic effect on the foliage 
profiles for the center scan of the red fir site 305, and the mean of five 
scans respectively. Although a good agreement (i.e., r = 0.977; 
RMSE = 0.0025 for center scan, and r = 0.991; RMSE = 0.0015 
for a mean of five scans) was found between the foliage profiles re- 
trieved before the topographic correction and after the topographic 
correction for these two cases, foliage profiles derived from the cen- 
ter scan and the mean foliage profiles derived from five scans dis- 
play an obvious topographic effect in terms of the shape of foliage 
profile and the topmost height. In particular, the peak at 10 m 
shrank greatly when the topographic correction was applied. In ad- 
dition, the shape shifted approximately 2-3 m, and resulted in a 
topmost height that was 2-3 m higher than foliage profiles without 
the local topographic correction. Fewer variations were observed 
for the mean foliage profiles from the five scans because the local 
topographic effects were somewhat averaged among the five EVI 
scans. For example, the peak at the 18 m level was significantly re- 
duced, and the topmost height also increased in the larger coverage. 

3.3. Foliage profile retrieved from EVI point cloud 

Fig. 8a, b demonstrates the effect of scanning perspective on foliage 
profiles derived from the EVI point cloud in two situations: before local 
topographic correction and after local topographic correction. The top- 
most height differences between these two situations increase over 
larger ranges, indicating the combined effect of the local topography 
and the spatial coverage. Another notable trend is that at closer ranges, 
the foliage profiles show smaller FAVD for the lower canopy (i.e. 
10-15 m), and larger FAVD for the upper canopy (30-40 m) except at 
20-m range. This can be explained by the fact that laser beams travel 
shorter distances before being intercepted at upper canopy level due 
to the high density of leaves, while the scattering events in the lower 
canopy tend to occur in the farther range due to lower density of leaves 
at this level. As a result, the horizontal ring setup tends to weight the 
upper canopy at closer ranges while for the lower canopy this effect is 
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Fig. 7. Topographic effect on foliage profiles derived from EVI full-waveform data, a) Center scan; b) mean of five scans. With slope, before slope correction. Without slope, after 
slope correction. 


dominated by scattering events occurring at farther range. The reason 
that this trend does not hold for 20 m range is probably due to a greatly 
clumped condition within the limited spatial coverage because EVI is 
usually located in a relatively open space. 

The topmost heights only vary very slightly when local topogra- 
phy is not considered. This is because the tree topmost heights are 
higher in the downslope regions of red fir site 305 and this tends to 
mask the local topographic effect somewhat. When the local topogra- 
phy is removed, the higher trees in the downslope are represented 
more appropriately. 

3.4. Foliage profile comparison 

Figs. 9 and 1 0 compare the foliage profiles retrieved from the follow- 
ing three sources: EVI full-waveform, EVI point cloud and LVIS. Two 



cases are considered for each data source: before a local topographic 
correction and after a local topographic correction. For most horizontal 
ranges, the foliage profiles retrieved from LVIS reveal higher topmost 
canopy heights than those retrieved from EVI. This reflects the fact 
that the airborne LVIS system, looking down from above, tends to char- 
acterize the canopy tips of the dominant trees more accurately while 
the terrestrial lidar system may miss the canopy tips of many of the 
dominant trees due to occlusion. This of course depends on the location 
of scanning point, and the forest tree density. In addition, LVIS-derived 
foliage profiles all demonstrate a higher upper canopy peak (at a height 
of 30-40 m) than those retrieved from EVI point cloud, due to the fact 
that LVIS recorded footprints illuminated within a d= 7° field of view 
while a nadir view was assumed in deriving foliage profile from EVI 
point cloud. Compared with the LVIS-derived foliage profiles, the up- 
ward looking EVI system samples the lower part of the forest stand 



Fig. 8. Foliage profiles by spatial coverage area radii using the EVI point cloud, a) Before slope correction; b) after slope correction. 


a) 30 m b) 40 m c) 50 m d) 60 m 
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Fig. 9. Foliage profile comparisons at varying horizontal ranges from the center before slope correction. 


a) 30 m b) 40 m c) 5 0 m d) 60 m 
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Fig. 10. Local topography corrected foliage profile comparisons at varying horizontal ranges from the center after slope correction. 
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more accurately, and thus may miss some of the top canopy part in the 
farther ranges, retrieving a lower FAVD for the upper canopy. 

Therefore, foliage profiles derived from the EVI full-waveform data 
all show higher FAVD for the lower canopy. The difference is caused pri- 
marily by the fact that the gap fractions are calculated for the maximum 
range, and as a result, more interceptions are contributing to the re- 
trieved foliage profiles, leading to a higher peak in the lower part of 
the canopy. As the horizontal ranges increase, the foliage profiles re- 
trieved from EVI full waveform correlate better, with r increasing from 
0.74 at 30 m to 0.92 at 100 m of horizontal distance. 

Foliage profiles from EVI point clouds agreed well with those from 
the LVIS system at all horizontal distances, but particularly in the farther 
range (e.g., r = 0.987 at 100 m horizontal distances). Note that the fo- 
liage profiles from the EVI point cloud are more similar to those re- 
trieved from LVIS, and are therefore actually apparent foliage profiles 
that are based on limited sampled zeniths. In contrast, the foliage pro- 
files from the EVI full-waveform data and those from the LVIS differ 
not only because of their different viewing directions, but also because 


of differences in the regions that are sampled. These comparisons sug- 
gest that the different methodologies used for foliage profile retrieval 
play an important role in the variation among foliage profiles retrieved 
from varied sources. 

After slope correction, the foliage profiles retrieved using the three 
methods display a similar pattern, as shown in Fig. 1 0. However, it ob- 
viously improves the correlation of foliage profiles from the EVI full 
waveform and the EVI point cloud, implying a critical need to remove 
local topographic effects before making foliage profile comparisons. 
The lower correlation between the profiles retrieved from LVIS and 
those retrieved from EVI point clouds also emphasizes the importance 
of correcting the local topography in the retrieval of foliage profiles at 
the shorter horizontal ranges. 

Fig. 1 1 a, b shows the effect of the scanning perspective on the foliage 
profile retrieval. Relatively lower correlations were observed for the 
smaller horizontal ranges (i.e., distance from center is less than 30 m). 
This may be due in part to the inherent errors in identifying the center 
points with the GPS devices. The local topographic effect also plays a 



b) After local topographic correction 



Fig. 11 . Correlation coefficient comparisons at varying horizontal ranges, a) Before local topography correction; b) after local topography correction. 


340 


E Zhao et al. / Remote Sensing of Environment 136 (2013) 330-341 


role. After the slope correction, the foliage profiles retrieved from the 
EVI full wave forms show a better correspondence with LVIS data than 
before the topographic correction, implying that the effects of topogra- 
phy at between-footprint scales had been removed, and consideration 
of the within-footprint topographic effects may provide more accurate 
foliage profiles from LVIS system. 

4. Discussion and conclusions 

This paper assesses the agreement between foliage profiles re- 
trieved from the airborne LVIS lidar scanning system with those re- 
trieved from the terrestrial Echidna lidar scanning system over 
varying horizontal distances from the site center, and over complex 
local topography. Airborne or spaceborne lidar scanning systems are 
attractive for mapping foliage profiles at regional and even continental 
scales because they can achieve larger spatial coverage, but they are 
limited by only capturing a near-nadir view and by an inability to 
characterize the detailed local topography. The comparison of foliage 
profiles from the lower spatial resolution, downward-looking air- 
borne LVIS lidar scanning system with the very highly sampled, 
upward-looking terrestrial EVI scanning system provides a better un- 
derstanding of the limitations of airborne lidar systems, as well as the 
mechanisms leading to the mismatch. 

This study reveals that view perspective, local topography, and spa- 
tial coverage are the main external factors that cause a mismatch of fo- 
liage profiles from these two lidar systems. We investigated the local 
topographic effect and spatial coverage effect on foliage profile re- 
trievals from three data sources: the EVI full-waveform data, the EVI 
point cloud, and the LVIS. We found that the correct characterization 
of topmost height was strongly affected by local topography, while 
the lower canopy and middle canopy correlated well regardless of 
whether the data had been topographically corrected or not. This find- 
ing suggests that height metrics retrieved from the LVIS system may 
need to consider the effect of local topography more carefully. With re- 
spect to the effect of scanning perspective, we have shown that in mak- 
ing direct comparisons of foliage profiles retrieved from small areas, the 
areas need to be closely matched, given GPS errors. With increasing area 
of spatial coverage (increasing horizontal range), the forests are sam- 
pled more evenly, and this results in profiles that are more similar. 

Our work shows the benefits of using the detailed description of forest 
canopy structure retrieved from terrestrial lidar systems in mapping fo- 
liage profiles over larger areas. As suggested by Ni-Meister et al. (2008), 
using a Geometric-Optical and Radiative-Transfer (GORT) model, relating 
both above-canopy and below-canopy lidar measurements may be 
expected to improve retrievals of forest canopy structure information, 
such as the foliage profile. With the help of terrestrial lidar, the negative 
impact of topography and the limited sampling region of foliage profiles 
can be minimized. 

Some similar findings have previously been reported by Hilker et al. 
(2010), who found good agreement between foliage profiles measured 
with EVI and airborne small-footprint foliage profiles in four plots of bo- 
real forest stands, although the airborne foliage profiles showed gener- 
ally less foliage in the middle and lower canopies. 

Our study has explored the comparison of foliage profiles between 
terrestrial and airborne large-footprint lidar systems by taking the effect 
of local topography, view perspective, and spatial coverage into account 
at a site dominated by sparsely distributed trees. 

The good agreement of foliage profile retrievals between EVI point 
cloud and the LVIS does not indicate that airborne measured foliage 
profiles are more accurate than those from EVI full waveform. Rather 
it suggests that the differences in the foliage profiles derived from 
the EVI full waveform data and LVIS system are mostly due to the 
way these two systems sample the forest. Although airborne lidar sys- 
tems have proven useful in retrieving foliage area profiles, they are 
limited by their nadir view; moreover, they must of necessity assume 
a predetermined leaf orientation (G function) or assume random leaf 


orientation, while the full-waveform EVI data provides a pathway to 
correct for G function effects. With a fine spatial resolution and a 
small beam size, terrestrial lidar systems complement the strengths 
of the airborne lidars by making a detailed characterization of the 
crowns from a small field site, thereby serving as a calibration tool 
and providing localized tuning information for future airborne and 
spaceborne lidar missions. 
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